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Abstract 

Anisotropic Kepler problem is investigated by perturbation method in both classi- 
cal and quantum mechanics. In classical mechanics, due to the singularity of the 
potential, global diffusion in phase space occurs at an arbitrarily small perturbation 
parameter. In quantum mechanics, the singularity induces a large transition am- 
plitude between quasi degenerate eigen states, which generically decays as Tl in the 
semi-classical limit. 
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The anisotropic Kelper problem (AKP) is a 2-dimensional Hamiltonian system defined by 

H = lipl+Pl)- , ^\ , (1) 

with fi > u > and ;Ui^ = 1 (For a detailed discussion of its dynamics, see ^ ). While the 
system reduces to the Kepler problem when /x = 1, the existence of chaos was rigorously proved 
by Gutzwiller and Devaney in the case of fi/v > 9/8 [||, ^] • The emergence of irregular dynamics 
dose not follow the conventional KAM scenario due to the singularity at the origin. In fact, the 
extremely unstable motion in the vicinity of the origin played a central role in Gutzwiller and 
Devaney's proof. The dynamical implication of the singularity when // — > 1 is, however, still not 
clear. When /x ~ 1, (1) can be rewritten as 

1/9 9n 1 cos 26 ^ 

H = -{vI+pI) ^ = ^0 + ei^, 2 

2 y r r 



where r = \J ^ and Q = ian~^{y/x) while e = (// — l)/2 appears as a natural perturbation 
parameter. In this letter we shall study the classical and quantum AKP by perturbation method. 
We hope this will provide a deeper understanding of the quantum-classic correspondence in non- 
KAM system. 

We begin with the classical mechanics. Kepler problem is a maximally integrable system. 
All its orbits are closed with period T = 2tt(—2E)~^/'^ , where E = Hq < 0. This global 
periodic motion will be destroyed by an infinitely small perturbation. The first order effect of 
a slight perturbation can be described by the drift of closed orbit, i.e., the slow variation of its 
parameters. Therefore, at this level of approximation, we can consider the motion in the space 
of all closed Kepler orbits (orbital space ). 

The orbital space can be conveniently constructed by taking advantage of the so(3) dynamical 
symmetry. Let 



Ji = (cos 6* - pyJs)/V-2E, J2 = {sm9 +pxJ3)/V-2E, and J3 = xpy - yp^. (3) 

(Ji, J2, J3) form a so(3) algebra, i.e., {Ji,Jj} = eijkJk,i,j = 1,2,3 and Hq = ^(fi+f^+fi) ■ 
There is a one-to-one correspondence between J' = (Ji, J2, J3) and close Kepler orbits. Specif- 



ically, write {Ji,J2-,Jz) = (sin a cos /3, sin a sin /3, cos a), the corresponding Kepler orbit in 
configuration space is defined as 

r = ^ (4) 

1 — sin a 008(6* — /3) 

while the sign of J3 determines its direction. Therefore, the orbital space is coordinated by 
iJi,J2,J3) (or iE,a,P)). 

The variation of Jj in one Kepler period is given by 

Ji{t + T) = Ji(t) + eT{Ji,Fo{J)} + o{e), i = 1,2,3, (5) 

where Fq is the time average of F over closed Kcplcr orbit, 

Fo = FoiJ) = - / F{r,9)—de = 2E—^ ^cos2/?. (6) 

r Jo Ksl l + |cosa| 

When e ^ 0, Eq. (5) can be approximated by a differential equation, 

^ = e{Ji,Fo{J)}, i = 1,2,3. (7) 

{Ho,Fq} = implies that the motion in the orbital space is confined within a sphere Se with 
Hq = E = const.. Restricting so(3) Poisson structure on Se induces a natural sympletic form 
U2 = -^^§d(3 A da, which is, up to a constant, the ordinary area element. 

As a two-dimensional Hamiltonian system, the dynamics on Se can be easily determined by 
the contour chart of the effective Hamiltonian Fq. There exist six fixed points on the sphere. 
The poles (a = 0, vr), which represent the two circular Kepler orbits, are unstable while the 
four on equator (a = it [2, (3 = 0, ib7r/2, tt), which corresponds to the linear orbits on the x or 
y axis, are stable. The remaining orbits are either the heteroclinic orbits (/3 = ±7r/4, ±37r/4) 
that connect the poles or periodic orbits surrounding one of the stable fixed points. In other 
words, besides separatrixs, the classical motion consists of four islands centered respectively at 
the degenerate Kepler orbits on the x and y axis. 

It should be pointed out that the simple picture given by perturbation analysis is not correct 
at the vicinity of a = 7r/2, i.e., the collision orbits, where eF is not bounded. In fact, the first 



order derivative of Fq is in general not continuous at equator. This non-smoothness manifests 
the non-perturbative nature of the motion near the origin. Noticing that all orbits on Se except 
for the poles cross a = 7r/2, we conclude that the global dynamics of AKP cannot be described 
by perturbation with respect to Kepler problem even at the limit e — 0. Our numerical study 
show that the collision orbits provide a passage-way for global diffusion in the orbital (and hence 
phase) space (Fig. 1). 

Now we turn to quantum mechanics. The so(3) symmetry of quantum planar Kepler problem 
is constructed in a way similar to its classical analog. Specifically, 

^1= [cos 6* - ^{pyJs + J3Pj/)]/\/-2i^o, 

< J2= [sme+l{p^J3 + Jsp^)]/V^2H^, (8) 

[Jj, Jfe] = iSjkiJi and Hq = ^^^"^1^ (h = 1). Let {|n, m >: m < \n\,n = 0, 1, ...} be the standard 
so(3) orthonormal set, i.e., J^\n,m >= m\n,m > and J'^\n,m >= n{n+l)\n,m >. \n,m > is an 
eigenstate of Hq with energy E = En = 2(^n+ly wave function in coordinate representation 
given by 

1 lin + ImlV 4sl"^l 1 p*™^ 

*»■•»<•■•'" = '-""^V in-Hi!(2„ + l)3/2 -'^-^'-" + ^1"! + ^'TS 

where s = 4r/ (2n -|- 1) and F(a, b, x) is the confluent hypergeometric function. 

The first order effect of a slight perturbation is the mixing of states with definite n. In the 
interaction representation, the long-time evolution when e — >^ is given by an effective Shordinger 
equation 

i-^qi = (10) 

UT 

where F is the restriction of F in eigen spaces of Hq, i.e., 

< n', m'\F\n, m >= Sn,n' < ^ Tn'|F|n, m > . (11) 
Eq. (10) can be regarded as the quantum counterpart of Eq. (6) and hence Fq is the classical 
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correspondence of F. The matrix elements of F is readily evaluated in coordinate representation, 



,,— , ^ in + mi)\{n — mi)\ . 

<n,m \F\n,m >= E„,J—— ^-fym' -m\,2, (12) 

y (n + m2)!(n — 777-2)! ' ' 

where ttii = min{|r77|, |r77'|} and 7772 = max{|r77|, |777'|}. 

By diagonalizing F in each 277 + 1-dimensional subspace, we can study the classical-quantum 
correspondence in the framework of perturbation theory. Notice that exp(7§ J3)F-|-Fexp(7^ J3) = 
0, the spectrum of F is symmetric with respect to 0. We shall focus on the positive part, which 
corresponds to the classical islands centered at fixed points (a,/?) = (7r/2,7r/2) and (7r/2, 37r/2). 

The classical orbits centered at fixed point (7r/2,7r/2) contribute to the spectrum of F ac- 
cording to the semi-classical quantization rule. 



J- [ co2 = k + Hs/4, A; = 0,1,... (13) 

where a{X) C Se is the region enclosed by orbit with Fo{E,a,P) = A and the Maslov index 
Us = 2. For convenience, we rescale E = — ^ and H = 2n+i so that |Fo| < 1 and Eq. (13) yields 

5(77, A,) ^ ^(cos-1 A, + ^^ii^) = k + \. (14) 

The orbits surrounding point (7r/2,37r/2) gives the same contribution so that the semi-classical 
spectrum consists of 2-fold degeneracies. Eq. (14) is a good approximation of the exact spectrum 
(Table 1). 

The degeneracy predicted by the semi-classical rule is not exact due to quantum tunneling 
between the two quasi static states each of them classically corresponds to one periodic orbit. 
The splitting of spectrum (AA) is connected with the transition amplitude accumulated in one 
Kepler period (^r) by At = eAXT/2?i. In conventional quantum system, such as a particle 
confined in double-well potential, AA decays as ~ ?P e-xp{—S/h) when h 0, which vanishes 
faster than any power of h. Fig. 2 shows AA^c = AA/^^ for 77 = 100 and 1000, from which we 
can see that A A ~ fi^ (or At ~ h) in the semi-classical limit. As a characteristic phenomenon 
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in in non-smooth system|Q], this power-law reflects the non-smoothness in Fq and hence the 
singularity in AKP. 

In summary, we have studied in this letter the Anisotropic Kepler problem by perturbation 
method. By using the dynamical symmetry of the unperturbed system, we obtained an effective 
Hamiltonian in both classical and quantum mechanics. We shown that the long-time evolution 
of most of the classical orbits is non-perturbative due to the singularity at origin and global 
diffusion in phase space occurs even at an arbitrarily small perturbation parameter. In quantum 
mechanics, we shown that the characteristic phenomenon attributed the singularity is the the 
power-law 7i-dependence of tunneling amplitude between states which classically correspond to 
separated periodic orbits. 
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Table 



k 


A 


5(n,A) 


A 


5(n,A) 





0.6883 


0.34 


0.6512 


0.42 


1 


0.3335 


1.52 


0.3264 


1.55 


2 


0.1944 


2.43 


0.1875 


2.48 


3 


0.0898 


3.50 


0.0856 


3.55 


4 


0.0307 


4.43 


0.0293 


4.46 



Tab.l Spectrum of F at ra = 10. 
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Figure 1: (a) Numerical calculation of one orbit at fi = 1.04. 5000 intersections are plotted on 
the Poincare surface of section defined at ^ = 0. with coordinate {R,P) = sign.{py){^/r, y/rpx), 
comparing with the regular structure implied by the first-order perturbation (b). 
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Figure 2: Scaled spectrum splitting (AAsc) for n = 100 and 1000. Except for few points near 
A = 0, (A, AAsc) is approximately located at two curves (determined by whether k is even or 
odd) irrespective of n. 
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